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Abstract 

O An analysis of I? 1 * 1 — >■ DK^ and -B 1 * 1 — >■ Dtt 1 * 1 decays is presented where the D 

meson is reconstructed in the four-body final state K^itTit + i{~ . Using LHCb data 
corresponding to an integrated luminosity of l.Ofb -1 , first observations are made of 
the suppressed ADS modes — > [k ± K t it + -k~]£)K ± and B^ — > [-K ± K T ir + Tr~]£)Tr ± 

^ with a significance of 5.1a and greater than lOu, respectively. Measurements of CP 

asymmetries and CP-conserving ratios of partial widths from this family of decays 
are also performed. The magnitude of the ratio between the suppressed and favoured 
B ± -> DK^ amplitudes is determined to be rf = 0.097 ± 0.011. 
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1 Introduction 



A measurement of the angle 7 (also denoted as 03 ) of the CKM Unitarity Triangle fl] 
in processes involving tree-level decays provides a Standard Model (SM) benchmark 
against which observables more sensitive to new physics contributions can be compared. 
Currently such comparisons are limited by the relatively large uncertainty (~ 12° [2]) on 
the determination of 7 in tree- level decays (3|[4]. More precise measurements are therefore 
required. 

A powerful strategy to measure 7 in tree-level processes is to study CP-violating 
observables in the decays B — > Dh , where D indicates a neutral charm meson which 
decays in a mode common to both D° and D° states, and h, the bachelor hadron, is 
either a kaon or a pion. In the case of B~ — > DK~, interference occurs between the 
suppressed b — > ucs and favoured b — > cus decay paths, and similarly for the charge 
conjugate decay. The magnitude of the interference is governed by three parameters: the 
weak-phase difference, 7, the CP-conserving strong-phase difference, 5^ , and the ratio of 
the magnitudes of the two amplitudes, r^. Similar interference effects occur in the case 
when the bachelor hadron is a pion, but additional Cabibbo suppression factors mean 
that the sensitivity to 7 is much diluted. Many possibilities exist for the D decay mode, 
including CP eigenstates j5) and self-conjugate three-body decays JH], which have both been 
exploited by LHCb in recent measurements (Tj[8j . Results of LHCb have also been presented 
making use of a similar strategy with B°/B° mesons |9|. Another option, termed the 



ADS' method in reference to its originators 10 , is to consider modes such as D — > K ± it T 
and to focus on the suppressed final-state B — > [tt ± K t ]dK ± , in which the favoured B ± 
decay is followed by a doubly Cabibbo-suppressed D decay, or the suppressed B ± decay 
precedes a favoured D decay. The amplitudes of such combinations are of similar total 
magnitude and hence large interference can occur, giving rise to significant CP-violating 
effects. In contrast, the interference in the favoured decay B ± — > [K ± ti t ]e ) K ± is low. 

In this Letter, a search is performed for the previously unobserved ADS decays 
B+- — > [7i ± K T 7i + 7i^]oh ± . The D decay is treated inclusively, with no attempt to separate 
out the intermediate resonances contributing to the four-body final state. LHCb has 
already presented an ADS study using B ± — > [tt^ K^ch^ decays, and many features of 
the current analysis are similar to those of the earlier paper (7j. In this study a total of 
seven observables is measured: the ratio of partial widths involving the favoured modes 



, K37T = T(B~ -» [K-7i + 7i-7i + } D K-) + T(B+ -» [K + 7i-7i + 71-} D K + ) 



two CP asymmetries, again involving the favoured modes 

T(B- -> [K-n+n+n-] D h~) - T(P+ -> [K+n-n+n~] D h+) 
h T(B- -> [K-7i+7i+7i-} D h-) + T(P+ -»■ [K+7i-7i+7i-] D h+) ' 1 ' 

and four charge-separated partial widths of the suppressed ADS mode relative to the 
favoured mode 

R K3.,± = r(g± -> [7c±KVn-] D h±) 



1 



The observables i?^ 37r,± carry the highest sensitivity to 7, r\ and 8\. They are related 



by the expression 10,11 



R K3n,± = (r | )2 + (r X3. )2 + 2r h r K^ K K^ + gKSn ± ^ (4) 

Here r^ 37r is the ratio of the magnitudes of the doubly Cabibbo-suppressed and Cabibbo- 
favoured D decay amplitudes, and 5^ 3n is the strong-phase difference between the ampli- 
tudes, averaged over the final-state phase space. The coherence factor K^ 37r accounts for 
possible dilution effects in the interference arising from the contribution of several interme- 
diate resonances in the D decay [ 1 1 1 . Information on these D decay parameters is available 
from external sources. Branching ratio measurements indicate that r^ 3 "" ~ 0.06 [12]. 
Studies of quantum-correlated DD pairs, performed at CLEO-c, yield 5^ 3n = (114 ^23)° 
and Kp 37T = 0.33 

-a23 The relatively low value of the coherence factor limits the 

sensitivity of R^ v,± to 7 and 5g, but does not hinder this observable in providing infor- 
mation on r\. Improved knowledge of r\ is valuable in providing a constraint which other 
—7- Dh^ analyses can benefit from. 



2 The LHCb detector and the analysis sample 

This analysis uses data, corresponding to an integrated luminosity of 1.0 fb _1 , collected 



by LHCb in 2011 at y\s = 7 TeV. The LHCb experiment 14 takes advantage of the 
high bb and cc production cross sections at the Large Hadron Collider to record large 
samples of heavy-hadron decays. It instruments the pseudorapidity range 2 < 77 < 5 of the 
proton-proton (pp) collisions with a dipole magnet and a tracking system that achieves 
a momentum resolution of 0.4 — 0.6% in the range 5 — 100 GeV/c. The dipole magnet 
can be operated in either polarity and this feature is used to reduce systematic effects 
due to detector asymmetries. In 2011, 58% of the data were taken with one polarity, 42% 
with the other. The pp collisions take place inside a silicon microstrip vertex detector that 
provides clear separation of secondary B ± vertices from the primary collision vertex (PV) 
as well as discrimination for tertiary D vertices. Two ring-imaging Cherenkov (RICH) 
detectors (15) with three radiators (aerogel, C 4 F 10 and CF 4 ) provide dedicated particle 
identification (PID), which is critical for the separation of B ± — >• DK ± and -^-D^ 
decays. 

A two-stage trigger is employed. First, a hardware-based decision is taken at a rate of 
up to 40 MHz. It accepts high transverse energy clusters in either an electromagnetic or 
hadron calorimeter, or a muon of high momentum transverse to the beam line (pr)- For 
this analysis, it is required that either one of the five tracks forming the B candidate 
points at a cluster in the hadron calorimeter, or that the hardware-trigger decision was 
taken independently of any of these tracks. A subsequent trigger level, implemented 
entirely in software, receives events at a rate of 1 MHz and retains ~ 0.3% of them. At 
least one track should have > 1-7 GeV/c and impact parameter (IP) x 2 with respect to 



The phase 5^ is given in the convention where CP\D°) = \D°). 
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the PV greater than 16. The IP y 2 is defined as the difference between the y 2 of the PV 
reconstructed with and without the considered track. In order to maximise efficiency at an 
acceptable trigger rate, a displaced vertex is selected with a decision tree algorithm that 
uses flight distance as well as fit quality, p^ and information on the IP with respect to the 
PV of the tracks. More information can be found in Ref. 16 . Full event reconstruction 



occurs offline, and a loose selection is run to reduce the size of the sample prior to final 
analysis. This selection consists of a decision tree algorithm similar to that used in the 
trigger, but in this case the entire decay chain is fully reconstructed and the selection 
benefits from the improved quality of the offline reconstruction. 

Approximately one million simulated signal events are used in the analysis as well 
as a sample of ~ 10 8 generic B q — > DX decays, where q G {u, d, s}. These samples are 
generated using Pythia 6.4 17 configured with parameters detailed in Ref. 18 . The 



EvtGen package 19 is used to generate hadronic decays, in which final state radiation is 
generated using the Photos package 20 . The interaction of the generated particles with 



the LHCb detector is simulated using the Geant4 toolkit 21 as described in Ref. 22 



3 Candidate selection and background rejection 

The reconstruction considers all B ± — >• Dh^ channels of interest. The reconstructed D 
candidate mass is required to be within ±25MeV/c 2 (~ 3.5a) of its nominal value [12]. 
The D daughter tracks are required to have p-r > 0.25GeV/c, while the bachelor track is 
required to satisfy 0.5 < pt < 10 GeV/c and 5 < p < 100 GeV/c. The tighter requirements 
on the bachelor track ensure that it resides within the kinematic coverage of the PID 
calibration samples acquired through the decay mode D* + — > D°tt + , D° — > K~tt + . Details 
of the PID calibration procedure are given in Sect. [4] Furthermore, a kinematic fit is 



performed to each decay chain 23 constraining both the B ± and D vertices to points in 
3D space, while simultaneously constraining the D candidate to its nominal mass. This 
fit results in a B ± mass resolution of 15MeV/c 2 , a 10% improvement with respect to the 
value prior to the fit. Candidates are retained that have an invariant mass in the interval 
5120 - 5750 MeV/c 2 . 

A boosted decision tree (BDT) discriminator |24|, implementing the GradientBoost 



algorithm 25 , is employed to achieve further background suppression. The BDT is trained 
using the simulated B ± — )■ Dh events together with a pure sample of combinatoric 
background candidates taken from a subset of the data in the invariant mass range 
5500 — 5800 MeV/ c 2 . The BDT considers a variety of properties associated with each signal 
candidate. These properties can be divided into two categories; (i) quantities common to 
both the tracks and to the D and B ± candidates, (ii) quantities associated with only the 
D and B ± candidates. Specifically, the properties considered in each case are as follows: 

(i) p, pt and IP y 2 ; 

(ii) decay time, flight distance from the PV, vertex quality, and the angle between the 
particle's momentum vector and a line connecting the PV to the particle's decay 
vertex. 
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In addition, the BDT also considers information from the rest of the event through an 
isolation variable that represents the imbalance of around the B ± candidate. The 
variable is defined as 

_ pr( B ) ~Ept 

where YIpt corresponds to the sum of over all tracks identified, within a cone of 
half-angle 1.5 in pseudorapidity and 1.5 rad in azimuthal angle, that are not associated 
with the signal B ± candidate. Since no PID information is used as input during the 
training, the BDT has very similar performance for both B ± — >• DK ± and B ± — > 
decay modes, with small differences arising from the variation in kinematics between the 
two. 

The optimal cut value on the BDT response is chosen by considering the combinatorial 
background level (b) in the invariant mass distribution of the favoured B — > Dtt^ 1 final 
state. The large signal peak in this sample is scaled by the predicted branching fraction 
of the suppressed mode for the case when the interference amongst the intermediate 
resonances of the D decay is maximally destructive. This provides a conservative estimate 
of the suppressed-sign signal yield (s). It is then possible to construct the quantity 
s/ys + b to serve as an optimisation metric. Assessment of this metric finds an optimal 
working point where a signal efficiency of ~ 85% is expected while rejecting > 99% of the 
combinatorial background. This same working point is used in selecting both suppressed 
and favoured final states. 

PID information is quantified as differences between the logarithm of likelihoods, InCh, 
under five mass hypotheses, h G {ir, K, p, e, /i} (DLL). The daughter kaon from the D 
meson decay is required to satisfy DLL^ = hivC^ — ln£„. > 2, while the daughter pions 
must have DLL^ < —2. A sample enriched in B ± — > DK ± decays is selected by requiring 
DLL^ > 4 for the bachelor hadron. Candidates failing this cut are retained in a separate 
sample, which is predominantly composed of B ± — > decays. 

Backgrounds from genuine B ± decays that do not involve a true D meson are suppressed 
by requiring the flight distance significance of the D candidate from the B^ vertex be greater 
than 2. The branching ratios of five-body charmless decays are currently unmeasured, and 
so the residual contamination from this source is estimated by assuming that the proportion 
of these decays passing the B — > [-K ± K T 7i + 7i^}i)h ± selection is the same as the proportion 
of three-body charmless decays passing the B ± —¥ K^dK^ selection reported in Ref. [7]. 
In the case of B — > n + n~ n + n~ K ± this assumption can be validated by removing the 
flight significance cut and inspecting the sideband above the D mass after adjusting the 
selection to isolate B ± — > [ii + 7i^7i + 7i^}£,K i: decays. A B — > n + n~ n + n~ K ± signal is 
observed with a magnitude compatible with that found when scaling the results of the 
analogous exercise performed with a B ± — > [h + 'k~]£)K ± selection, and is eliminated when 
the flight- significance cut is reinstated. Following these studies the residual contamination 
from five-body charmless decays is determined to be 2.2 ± 1.1 candidates for the suppressed 
B ± — > DK ± selection, and negligible for all other samples. 

Contamination involving misidentified charmonium decays is eliminated by considering 
the possible neutral combinations of the bachelor track and any one of the D daughter 
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tracks under the hypothesis that both tracks are muons. For those combinations where 
both tracks satisfy a loose muon PID requirement, the parent B ± candidate is vetoed if 
the invariant mass of this combination is within ±22MeV/c 2 of either the J/if) or ip(2S) 
mass 
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The suppressed signal sample suffers a potentially large cross-feed from favoured 
signal decays in which a K and 7r T from the D decay are misidentified as 7r =l= and K T , 
respectively. This contamination is reduced by vetoing any suppressed candidate whose 
reconstructed D mass, under the exchange of mass hypotheses between the daughter kaon 
and either of the two same-sign daughter pions, lies within ±15MeV/c 2 of the nominal 
D mass. For the measurement of R^ 3n,± this veto is also applied to the favoured mode. 
Study of the cross-feed contamination in the mass sidebands of the D candidates allows 
the estimate of the residual contamination in the signal region to be checked. The residual 
cross-feed after all selection requirements is estimated to be (7.1 ±3.1) x 10 -5 . 

The mass window for the D candidates is sufficiently tight to eliminate background 
arising from single-track misidentifications in the four-body decays D — > K + K + ti + tt~ and 
D — > 7T+ 7i~ 7i + 7i~ . The good performance of the RICH system ensures that the residual 
background from D — > K + K~7t + tt~ decays in which three tracks are misidentified is 
negligible. The contamination in the suppressed B — > [7r ± K T 7i + 7i^]£,K i: sample from 

—7- [K^ K®h~]dK ± with — > 7r + 7r~ is estimated from simulation to occur at the rate 
of (6.1 ± 1.7) x 1(T 5 . 

Only one candidate per event is retained for analysis. In the 0.8% of events that 
contain more than one candidate a choice is made by selecting the candidate with the 
best-quality B ± vertex. 

Using simulation it is found that the selection leads to an acceptance across the four- 
body phase space of the D decay that is uniform to a good approximation. This property 
is important as it means that the values of the coherence factor and strong-phase difference 



measured in Ref. 13 , which are integrated over all phase space, can be applied when 
interpreting the results of the current analysis. It is verified using simulation that the 
small non-uniformities that exist in the acceptance induce negligible bias in the effective 
value of these parameters. 



4 Signal yields and systematic uncertainties 

The observables of interest are determined with a binned maximum-likelihood fit to the 
invariant mass distributions of the selected B ± candidates. Distinguishing between B + 
and B~ candidates, favoured and suppressed decay topologies, and those that pass or 
fail the bachelor PID requirement imposed to select B ± — > DK^ decays, yields eight 
disjoint subsamples, which are fitted simultaneously. The total probability density function 
(PDF) used in the fit is built from four main sources representing the various categories of 
candidates in each subsample. 
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1. B ± Dn* 

In the subsamples failing the bachelor PID cut, a modified Gaussian function, 

f(m) oc exp ( (™ — ^- \ (6) 

\2a 2 + {m- H) 2 a LjR J 

describes the asymmetric function of peak of value fi and width a where a^(m < jj) 
and «ij(m > /i) parameterise the tails. True B ± — » /^vr 1 * 1 candidates that pass the 
PID cut are reconstructed as B — > DK ± . As these candidates have an incorrect 
mass assignment they form a displaced mass peak with a tail that extends to higher 
invariant mass. These candidates are modelled by the sum of two Gaussian PDFs, 
also altered to include tail components as in Eq. [6} All parameters are allowed to 
vary except the lower-mass tail which is fixed to the value found in simulation to 
ensure fit stability, and later considered amongst the systematic uncertainties. These 
shapes are considered identical for B + and B~ decays. 

2. B ± DK 

In the subsamples that pass the PID cut on the bachelor, the same modified Gaussian 
function as quoted in Eq. [6] is used. The peak value and the two tail parameters are 
identical to those of the higher B ± — > peak. The width is 0.95 ± 0.02 times 
the B ± — > Dti^ 1 width, as determined by a separate study of the favoured mode. 
Candidates failing the PID cut are described by a fixed shape that is obtained from 
simulation and later varied to assess the systematic uncertainty. 

3. Partially reconstructed b-hadron decays 

Partially reconstructed decays populate the invariant mass region below the B mass. 
Such candidates may enter the signal region, especially where Cabibbo-favoured 
B ± -> XDvr ± modes are misidentified as B ± — > DK ± . The large simulated sample 
of inclusive B q — > DX decays is used to model this background. After applying 



the selection, two non-parametric PDFs 26 are defined (for the and DK ± 



selections) and used in the signal fit for both the favoured and suppressed mode 
subsamples. 

In addition, partially reconstructed B® — > DK~tt + and A® — > [pK~7T + 7T~7T + ] A +h~ 
decays and their charge- conjugated modes are considered as background sources 
specific to the suppressed B ± — > DK ± and favoured mode subsamples, respectively. 
PDFs for both these sources of background are determined from simulation and 
smeared to match the resolution observed in data. 

The yield of these background components in each subsample varies independently 
in the fit, making no assumption of CP symmetry. 

4. Combinatoric background 

A linear approximation is adequate to describe the distribution across the invariant 
mass spectrum considered. A common shape parameter is used in all subsamples, 
though yields vary independently. 
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Table 1: Favoured {[Kttt^tt]^) and suppressed ([TrKTnr]nh) signal yields together with their 
corresponding statistical uncertainties. 



Mode 


B~ 


B+ 




20, 791 ± 232 


21, 054 ±235 


[K-KirTrjoK 


1,567± 57 


1,660± 60 




87 ± 11 


68 ± 10 




11 db 5 


29 ± 7 



The proportion of B ± — > Dh^ 1 passing or failing the PID requirement is determined from 
an analysis of approximately 20 million D decays reconstructed as D*^ — > Dir^, D — > 
K T 7i ± . The reconstruction is performed using only kinematic variables, and provides a 
high purity calibration sample of K and tt tracks which is unbiased for studies exploiting 
the RICH and is therefore made use of to measure the PID efficiency as a function of track 
momentum, pseudorapidity and number of tracks in the detector. Through reweighting the 
calibration spectra in these variables to match that of the candidates in the B — > 
peak, the effective PID efficiency of the signal is determined. This data-driven approach 
finds a retention rate on the bachelor track of 86.1% and 3.7% for kaons and pions, 
respectively. An absolute 1.0% systematic uncertainty on the kaon efficiency is estimated 
from simulation. The B ± — > Dn^ fit to data becomes significantly incorrect when the 
PID efficiency is varied outside the absolute range of ±0.2%, and so this value is taken as 
the systematic uncertainty for pions. 

Detection and production asymmetries are accounted for using the same procedure 
followed in Ref. [7], based on the measurement of the observed raw asymmetry of B ± — > 



J/ipK ± decays in the LHCb detector 27 . A detection asymmetry of (—0.5 ± 0.7)% is 
assigned for each unit of strangeness in the final state to account for the different interaction 
lengths of K~ and K + mesons. The equivalent asymmetry for pions is expected to be 
much smaller and (0.0 ±0.7)% is assigned. This uncertainty also accounts for the residual 
physical asymmetry between the left and right sides of the detector after summing both 
magnet-polarity data sets. Simulation of 6-hadron production in pp collisions suggests 
a small excess of B + over B~ mesons. A production asymmetry of (—0.8 ± 0.7)% is 
assumed in the fit such that the combination of these estimates aligns with the observed 



raw asymmetry of B — > J/ipK decays 27 



The signal yields for the favoured and suppressed B — >■ Dh^ 1 decays, after summing the 
events that pass and fail the bachelor PID cut, are shown in Table [T] Their corresponding 
invariant mass spectra, separated by the charge of the B candidate, are shown in Figs. [T] 
and [2] Plots of the combined B + and B~ suppressed-mode mass spectra are shown in 
Fig-i 

The R h 7r ' observables are related to the ratio of event yields by the relative efficiency, 
independent of PID effects, with which B ± — > DK ± and B^ — > Dtc^ 1 decays are recon- 
structed. This ratio is found to be 1.029 from a simulation study. A 2.4% systematic 
uncertainty, based on the finite size of the simulated sample, accounts for the imperfect 
modelling of the relative pion and kaon absorption in the tracking material. 



7 




5200 



5400 
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Figure 1: Invariant mass distributions of selected — > [K^Tr^TT^Tr^jnh^ 1 candidates, separated 
by charge. The left plots are B~ candidates, B + are on the right. In the top plots, the bachelor 
track passes the PID cut and the B^ candidates are reconstructed assigning this track the kaon 
mass. The remaining candidates are placed in the sample displayed on the bottom row and are 
reconstructed with a pion mass hypothesis. The dark (red) and light (green) curves represent the 
fitted B^ — > DK^ and B^ — > Dir^ components, respectively. The shaded contribution indicates 
partially reconstructed decays and the total PDF includes the combinatorial component. 



The fit is constructed such that the observables of interest are free parameters. To 
estimate the systematic uncertainties arising from the imperfect knowledge of several of 
the external parameters discussed above, the fit is performed many times varying each 
input by its assigned error. The resulting spread (RMS) in the value of each observable 
is taken as the systematic uncertainty on that quantity and is summarised in Table |2j 
Correlations between the uncertainties are considered negligible, so the total systematic 
uncertainty is the sum in quadrature of the individual components. 



S 



> 

c 



> 




5200 



5400 



5600 



5800 



5200 



5400 



5600 5800 

m(D/i ± ) [MeV/c 2 ] 



Figure 2: Invariant mass distributions of selected — > [ir^ K^^-k^]^^ decays, separated 
by charge. See the caption of Fig. [I] for a full description. The dashed line here represents the 
partially reconstructed, but Cabibbo favoured, B® — > DK~7r + , and charge-conjugated, decays 
where the pion is not reconstructed. The favoured mode cross- feed is included in the fit, but is 
too small to be seen. 



Table 2: Systematic uncertainties on the observables. Bachelor PID refers to the fixed efficiency 
for the bachelor track DLL^ requirement determined using the D*^ calibration sample. PDFs 
refers to the variations of the fixed shapes in the mass fit. Simulation refers to the use of 
simulation to estimate relative efficiencies of the signal modes, and also includes the contribution 
from the uncertainty in the residual background from charmless B decays. ^4instr. quantifies the 
uncertainty on the production, interaction and detection asymmetries. 



[xi(r 3 ] 


T)K3lT 
a K/7T 


^K3tt 








u K3tt,- 




Bachelor PID 


1.7 


0.2 


0.6 


0.4 


0.4 


0.02 


0.04 


PDFs 


1.2 


1.3 


4.4 


0.7 


0.9 


0.09 


0.08 


Simulation 


1.5 


0.1 


0.3 


0.1 


0.2 


0.01 


0.02 


-^■instr. 


0.0 


9.9 


17.1 


0.1 


0.1 


0.06 


0.06 


Total 


2.6 


10.0 


17.7 


0.8 


1.0 


0.11 


0.11 
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Figure 3: Charge-integrated invariant mass distributions of those candidates shown in Fig. [2] 
passing the — > [n^K^ir^TT^jDh^ 1 selection. 



5 Results and interpretation 

The results of the fit with their statistical and systematic uncertainties are 

= 0.077i ± 0.0017 ± 0.0026 



dK3it 
n K/n 

J[K3tt 



R 
R 



K3ir,- 
K 

K3-w,+ 
K 

R K3«,+ 



-0.029 
-0.006 
0.0072 
0.0175 



± 0.020 ± 0.018 
± 0.005 ± 0.010 



+ 0.0036 

- 0.0032 

+ 0.0043 

- 0.0039 



± 0.0008 
± 0.0010 



0.00417 t aoooso ± 0.00011 
0.00321 + ± 0.00011 



From these measurements, the quantities -Rads(/i) an( ^ ^ads(/i) can ^ e deduced. These 
are, respectively, the ratio of the suppressed to the favoured partial widths for the decays 
5 ± — > Dh^, averaged over the two charges, and the CP asymmetry of the suppressed 



10 



decay mode 



jjK3tt,— i jjK3tt, 
R ADS(K) = — ~ " = ± {]m27 



Afll {K) = ^kz ^^ = -0.42 ± 0.22 





2 


r>K3w - 
K 


- R K 


K K 


+ R K 


RK3TT- 


+ R^ + 




2 


R K3n ~ 


tdK3it,+ 

- fX 7r 




, tjK3-k,+ 



R ABS(n) = — 5 = °- 0037 ± °- 0004 

Af^ (7r) = ^ = 0.13 ± 0.10 . 

The displayed uncertainty is the combination of statistical and systematic contributions. 
Correlations between systematic uncertainties are taken into account in the combination. 
It can be seen that the observable ^ADS(if)' which is expected to manifest significant CP 
violation, differs from the CP-conserving hypothesis by around 2a. 

A likelihood ratio test is employed to assess the significance of the suppressed ADS 

signal yields. This has been performed calculating the quantity «/— 2 In 7^-, where £ s +b 

y ^s+b 

and £b ar e the maximum values of the likelihoods in the case of a signal-plus-background 
and background-only hypothesis, respectively. Significances of 5.7cr and greater than 
10<7 are determined for the modes B — > [TT ± K ± 7r + 7r~]r)K ± and B ± — > {'k ± K t 7t + 7t~}dtt ± , 
respectively. The former significance is found to reduce to 5.1a when the systematic 
uncertainties are included. 

The measured observables are used to infer a confidence interval for the value of the 
suppressed-to-favoured B — >■ DK amplitude ratio, rjj. The most probable value of 
is identified as that which minimises the y 2 calculated from the measured observables anc 
their predictions for the given value of r^. The prediction for R^ 3n,± is given by Eq. [4 
and similar relations exist for the other observables. Amongst the other parameters that 
determine the predicted values, r^, 8^ and 7 vary freely, but all the parameters of the D 
decay, notably the coherence factor and strong- phase difference, are constrained by the 
results in Ref. 13 . Subsequently, the evolution of the minimum \ 2 is inspected across 



the range (0.0 <rg < 0.2) and the difference A% 2 with respect to the global minimum is 
calculated. The probabilistic interpretation of the Ax 2 at each value of is evaluated by 
generating and fitting a large number (10 7 ) of pseudo-datasets around the local minimum. 
The variation of the pseudo-datasets is derived from the covariance matrix of the principal 
result. At a given fixed point in the range, a, with a value of A\ 2 a above the global 
minimum, the probability of obtaining the observed result is defined as the number of 
pseudo-experiments with A% 2 > A% 2 . By this frequentist technique it is found that the 
result for r§ has a non-Gaussian uncertainty, so the "lcr" and "2a" intervals, respectively, 
are given as 

r f = 0.097 ±0.011 [68.3% CL] and t° f 2 l [95.5% CL] . 
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The measurements do not allow significant constraints to be set on the other underlying 
physics parameters. 



6 Conclusions 

A search has been performed for the ADS suppressed modes B ± — > [k ± K t 'k + 'k~]e>K ± 
and B ± [^K^ti+tt^d^ using l.Ofb -1 of data collected by LHCb in 2011. First 
observations have been made of both decays, with a significance of 5.1a and greater than 
10cr, respectively. Measurements have been made of the observables R^%, A^ 3lT and 

i?^ 37r,± , as well as the derived parameters R^BS(h) anc ^ ^ads(/i) » which relate the partial 
widths of the B ± — > Dh^ 1 (h = K, tc) family of decays. From these observables it is deduced 
that = 0.097 ± 0.011, where is the ratio of the absolute values of the suppressed and 
favoured B ± — > DK ± amplitudes. These results will improve knowledge of the Unitarity 
Triangle angle 7 when they are combined with other B ± — > DK ± measurements exploiting 
different D decay modes. 
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